Germinal Center Dark Zone harbors ATR-dependent determinants of T-cell exclusion that are also identified in aggressive lymphoma.

The germinal center (GC) dark zone (DZ) and light zone (LZ) regions spatially separate expansion and diversification from selection of antigen-specific B-cells to ensure antibody affinity maturation and B cell memory. The DZ and LZ differ significantly in their immune composition despite the lack of a physical barrier, yet the determinants of this polarization are poorly understood. This study provides novel insights into signals controlling asymmetric T-cell distribution between DZ and LZ regions. We identify spatially-resolved DNA damage response and chromatin compaction molecular features that underlie DZ T-cell exclusion. The DZ spatial transcriptional signature linked to T-cell immune evasion clustered aggressive Diffuse Large B-cell Lymphomas (DLBCL) for differential T cell infiltration. We reveal the dependence of the DZ transcriptional core signature on the ATR kinase and dissect its role in restraining inflammatory responses contributing to establishing an immune-repulsive imprint in DLBCL. These insights may guide ATR-focused treatment strategies bolstering immunotherapy in tumors marked by DZ transcriptional and chromatin-associated features.

Spatially-resolved transcriptomics reveal macrophage heterogeneity and prognostic significance in diffuse large B-cell lymphoma.

Macrophages are abundant immune cells in the microenvironment of diffuse large B-cell lymphoma (DLBCL). Macrophage estimation by immunohistochemistry shows varying prognostic significance across studies in DLBCL, and does not provide a comprehensive analysis of macrophage subtypes. Here, using digital spatial profiling with whole transcriptome analysis of CD68+ cells, we characterize macrophages in distinct spatial niches of reactive lymphoid tissues (RLTs) and DLBCL. We reveal transcriptomic differences between macrophages within RLTs (light zone /dark zone, germinal center/ interfollicular), and between disease states (RLTs/ DLBCL), which we then use to generate six spatially-derived macrophage signatures (MacroSigs). We proceed to interrogate these MacroSigs in macrophage and DLBCL single-cell RNA-sequencing datasets, and in gene-expression data from multiple DLBCL cohorts. We show that specific MacroSigs are associated with cell-of-origin subtypes and overall survival in DLBCL. This study provides a spatially-resolved whole-transcriptome atlas of macrophages in reactive and malignant lymphoid tissues, showing biological and clinical significance.

ANGPTL3 Downregulation Increases Intracellular Lipids by Reducing Energy Utilization.

BACKGROUND: ANGPTL3 (angiopoietin-like protein 3) is a circulating protein with a key role in maintaining lipoprotein homeostasis. A monoclonal antibody against ANGPTL3 is an approved and well-tolerated treatment to reduce lipoproteins in familial hypercholesterolemia homozygotes. However, the reduction of hepatic ANGPTL3 synthesis using an antisense oligonucleotide unexpectedly resulted in a dose-dependent increase in liver lipid content and circulating transaminases, resulting in the termination of the clinical trial. Meanwhile, the use of silencing RNAs remains an area of active investigation. Our study sought to investigate whether intracellular downregulation of ANGPTL3 may lead to a primary increase in neutral lipids within the hepatocyte. METHODS: We downregulated ANGPTL3 by silencing RNA in primary human hepatocytes 3-dimensional spheroids, HepG2/LX-2 3-dimensional spheroids, and in HepG2, Hep3B2, and Huh7 cultured in 2 dimensions. RESULTS: ANGPTL3 downregulation increased neutral lipids in all models investigated. Interestingly, ANGPTL3 induced lower intracellular deiodinase type 1 protein levels resulting in a reduction in beta-oxidation and causing an increase in triglycerides stored in lipid droplets. CONCLUSIONS: In conclusion, intracellular ANGPTL3 downregulation by silencing RNA led to an increase in triglycerides content due to a reduction in energy substrate utilization resembling a primary intracellular hepatocyte hypothyroidism.

No Clear Clustering Dysbiosis from Salivary Microbiota Analysis by Long Sequencing Reads in Patients Affected by Oral Squamous Cell Carcinoma: A Single Center Study.

BACKGROUND: Advancements in DNA sequencing technology have facilitated the assessment of the connection between the oral microbiome and various diseases. The aim of the present study was to investigate the salivary microbiota composition employing for the first time in the literature the Oxford Nanopore Technology in patients affected by oral squamous cell carcinoma (OSCC). METHODS: Unstimulated saliva samples of 31 patients were collected (24 OSCC patients and 7 controls). DNA was extracted using the QIAamp DNA Blood Kit and metagenomic long sequencing reads were performed using the MinION device. RESULTS: In the OSCC group, 13 were males and 11 were females, with a mean age of 65.5 ± 13.9 years; in the control group, 5 were males and 2 were females, with a mean age of 51.4 ± 19.2 years. The border of the tongue was the most affected OSCC site. The microorganisms predominantly detected in OSCC patients were Prevotella, Chlamydia, Tissierellia, Calothrix, Leotiomycetes, Firmicutes and Zetaproteobacteria. CONCLUSIONS: This study confirmed the predominance of periodontopathic bacteria in the salivary microbiome in the OSCC group. If a direct correlation between oral dysbiosis and OSCC onset was proven, it could lead to new prevention strategies and early diagnostic tools.

Pseudotemporal ordering of spatial lymphoid tissue microenvironment profiles trails Unclassified DLBCL at the periphery of the follicle.

We have established a pseudotemporal ordering for the transcriptional signatures of distinct microregions within reactive lymphoid tissues, namely germinal center dark zones (DZ), germinal center light zones (LZ), and peri-follicular areas (Peri). By utilizing this pseudotime trajectory derived from the functional microenvironments of DZ, LZ, and Peri, we have ordered the transcriptomes of Diffuse Large B-cell Lymphoma cases. The apex of the resulting pseudotemporal trajectory, which is characterized by enrichment of molecular programs fronted by TNFR signaling and inhibitory immune checkpoint overexpression, intercepts a discrete peri-follicular biology. This observation is associated with DLBCL cases that are enriched in the Unclassified/type-3 COO category, raising questions about the potential extra-GC microenvironment imprint of this peculiar group of cases. This report offers a thought-provoking perspective on the relationship between transcriptional profiling of functional lymphoid tissue microenvironments and the evolving concept of the cell of origin in Diffuse Large B-cell Lymphomas.

Programmed cell death 1 genetic variant and liver damage in nonalcoholic fatty liver disease.

BACKGROUND AND AIMS: Programmed cell death 1/programmed cell death-ligand 1 (PD-1/PDL-1) axis has been reported to modulate liver inflammation and progression to hepatocellular carcinoma (HCC) in patients with nonalcoholic fatty liver disease (NAFLD). Here, we examined whether the PDCD1 variation is associated with NAFLD severity in individuals with liver biopsy. METHODS: We examined the impact of PDCD1 gene variants on HCC, as robust severe liver disease phenotype in UK Biobank participants. The strongest genetic association with the rs13023138 G>C variation was subsequently tested for association with liver damage in 2889 individuals who underwent liver biopsy for suspected nonalcoholic steatohepatitis (NASH). Hepatic transcriptome was examined by RNA-Seq in a subset of NAFLD individuals (n = 121). Transcriptomic and deconvolution analyses were performed to identify biological pathways modulated by the risk allele. RESULTS: The rs13023138 C>G showed the most robust association with HCC in UK Biobank (p = 5.28E-4, OR = 1.32, 95% CI [1.1, 1.5]). In the liver biopsy cohort, rs13023138 G allele was independently associated with severe steatosis (OR 1.17, 95% CI 1.02-1.34; p = .01), NASH (OR 1.22, 95% CI 1.09-1.37; p < .001) and advanced fibrosis (OR 1.26, 95% CI 1.06-1.50; p = .007). At deconvolution analysis, rs13023138 G>C allele was linked to higher hepatic representation of M1 macrophages, paralleled by upregulation of pathways related to inflammation and higher expression of CXCR6. CONCLUSIONS: The PDCD1 rs13023138 G allele was associated with HCC development in the general population and with liver disease severity in patients at high risk of NASH.

Salivary Microbiota Composition in Patients with Oral Squamous Cell Carcinoma: A Systematic Review.

BACKGROUND: Oral squamous cell carcinoma (OSCC) is one of the most prevalent cancers worldwide. Despite recent advances in diagnosis and treatment, in recent years, an increase in the incidence of OSCC has been registered, and the mortality rate is still high. This systematic review aims to identify a potential association between the composition of salivary microbiota and OSCC. MATERIALS AND METHODS: The protocol for this study was designed following the PRISMA guidelines. Records were identified using different search engines (e.g., Medline/PubMed). Observational studies, in human subjects with histological diagnosis of OSCC, concerning the analysis of salivary microbiota, were selected. RESULTS: Eleven papers were included. The salivary microbiomes of 1335 patients were analysed (n.687 OSCC and n.648 controls). Due to the great heterogeneity of the studies, it was not possible to profile a specific microbiota associated with OSCC. However, periodontal pathogens were the most common bacteria detected in patients with OSCC (i.e., Fusobacterium, Prevotella). CONCLUSIONS: Although there are evident alterations in the salivary microbiota composition in OSCC patients, it is still a challenge to identify a specific microbiota pattern in OSCC patients. If the associations between specific salivary microorganisms and OSCC are confirmed, microbiome analysis could be a useful tool for the screening and follow-up of patients affected by OSCC.

Spatial transcriptome of a germinal center plasmablastic burst hints at MYD88/CD79B mutants-enriched diffuse large B-cell lymphomas.

The GC reaction results in the selection of B cells acquiring effector Ig secreting ability by progressing toward plasmablastic differentiation. This transition is associated with exclusion from the GC microenvironment. The aberrant expansion of plasmablastic elements within the GC fringes configures an atypical condition, the biological characteristics of which have not been defined yet. We investigated the in situ immunophenotypical and transcriptional characteristics of a nonclonal germinotropic expansion of plasmablastic elements (GEx) occurring in the tonsil of a young patient. Compared to neighboring GC and perifollicular regions, the GEx showed a distinctive signature featuring key regulators of plasmacytic differentiation, cytokine signaling, and cell metabolism. The GEx signature was tested in the setting of diffuse large B-cell lymphoma (DLBCL) as a prototypical model of lymphomagenesis encompassing transformation at different stages of GC and post-GC functional differentiation. The signature outlined DLBCL clusters with different immune microenvironment composition and enrichment in genetic subtypes. This report represents the first insight into the transcriptional features of a germinotropic plasmablastic burst, shedding light into the molecular hallmarks of B cells undergoing plasmablastic differentiation and aberrant expansion within the noncanonical setting of the GC microenvironment.

A multivariate statistical test for differential expression analysis.

Statistical tests of differential expression usually suffer from two problems. Firstly, their statistical power is often limited when applied to small and skewed data sets. Secondly, gene expression data are usually discretized by applying arbitrary criteria to limit the number of false positives. In this work, a new statistical test obtained from a convolution of multivariate hypergeometric distributions, the Hy-test, is proposed to address these issues. Hy-test has been carried out on transcriptomic data from breast and kidney cancer tissues, and it has been compared with other differential expression analysis methods. Hy-test allows implicit discretization of the expression profiles and is more selective in retrieving both differential expressed genes and terms of Gene Ontology. Hy-test can be adopted together with other tests to retrieve information that would remain hidden otherwise, e.g., terms of (1) cell cycle deregulation for breast cancer and (2) "programmed cell death" for kidney cancer.

RNA Sequencing of Primary Cutaneous and Breast-Implant Associated Anaplastic Large Cell Lymphomas Reveals Infrequent Fusion Transcripts and Upregulation of PI3K/AKT Signaling via Neurotrophin Pathway Genes.

Cutaneous and breast implant-associated anaplastic large-cell lymphomas (cALCLs and BI-ALCLs) are two localized forms of peripheral T-cell lymphomas (PTCLs) that are recognized as distinct entities within the family of ALCL. JAK-STAT signaling is a common feature of all ALCL subtypes, whereas DUSP22/IRF4, TP63 and TYK gene rearrangements have been reported in a proportion of ALK-negative sALCLs and cALCLs. Both cALCLs and BI-ALCLs differ in their gene expression profiles compared to PTCLs; however, a direct comparison of the genomic alterations and transcriptomes of these two entities is lacking. By performing RNA sequencing of 1385 genes (TruSight RNA Pan-Cancer, Illumina) in 12 cALCLs, 10 BI-ALCLs and two anaplastic lymphoma kinase (ALK)-positive sALCLs, we identified the previously reported TYK2-NPM1 fusion in 1 cALCL (1/12, 8%), and four new intrachromosomal gene fusions in 2 BI-ALCLs (2/10, 20%) involving genes on chromosome 1 (EPS15-GNG12 and ARNT-GOLPH3L) and on chromosome 17 (MYO18A-GIT1 and NF1-GOSR1). One of the two BI-ALCL samples showed a complex karyotype, raising the possibility that genomic instability may be responsible for intra-chromosomal fusions in BI-ALCL. Moreover, transcriptional analysis revealed similar upregulation of the PI3K/Akt pathway, associated with enrichment in the expression of neurotrophin signaling genes, which was more conspicuous in BI-ALCL, as well as differences, i.e., over-expression of genes involved in the RNA polymerase II transcription program in BI-ALCL and of the RNA splicing/processing program in cALCL.

Constitutive PSGL-1 Correlates with CD30 and TCR Pathways and Represents a Potential Target for Immunotherapy in Anaplastic Large T-Cell Lymphoma.

Due to the high expression of P-selectin glycoprotein ligand-1 (PSGL-1) in lymphoproliferative disorders and in multiple myeloma, it has been considered as a potential target for humoral immunotherapy, as well as an immune checkpoint inhibitor in T-cells. By investigating the expression of SELPLG in 678 T- and B-cell samples by gene expression profiling (GEP), further supported by tissue microarray and immunohistochemical analysis, we identified anaplastic large T-cell lymphoma (ALCL) as constitutively expressing SELPLG at high levels. Moreover, GEP analysis in CD30+ ALCLs highlighted a positive correlation of SELPLG with TNFRSF8 (CD30-coding gene) and T-cell receptor (TCR)-signaling genes (LCK, LAT, SYK and JUN), suggesting that the common dysregulation of TCR expression in ALCLs may be bypassed by the involvement of PSGL-1 in T-cell activation and survival. Finally, we evaluated the effects elicited by in vitro treatment with two anti-PSGL-1 antibodies (KPL-1 and TB5) on the activation of the complement system and induction of apoptosis in human ALCL cell lines. In conclusion, our data demonstrated that PSGL-1 is specifically enriched in ALCLs, altering cell motility and viability due to its involvement in CD30 and TCR signaling, and it might be considered as a promising candidate for novel immunotherapeutic approaches in ALCLs.

miR-1207-5p Can Contribute to Dysregulation of Inflammatory Response in COVID-19 via Targeting SARS-CoV-2 RNA.

The present study focuses on the role of human miRNAs in SARS-CoV-2 infection. An extensive analysis of human miRNA binding sites on the viral genome led to the identification of miR-1207-5p as potential regulator of the viral Spike protein. It is known that exogenous RNA can compete for miRNA targets of endogenous mRNAs leading to their overexpression. Our results suggest that SARS-CoV-2 virus can act as an exogenous competing RNA, facilitating the over-expression of its endogenous targets. Transcriptomic analysis of human alveolar and bronchial epithelial cells confirmed that the CSF1 gene, a known target of miR-1207-5p, is over-expressed following SARS-CoV-2 infection. CSF1 enhances macrophage recruitment and activation and its overexpression may contribute to the acute inflammatory response observed in severe COVID-19. In summary, our results indicate that dysregulation of miR-1207-5p-target genes during SARS-CoV-2 infection may contribute to uncontrolled inflammation in most severe COVID-19 cases.

An improvement of ComiR algorithm for microRNA target prediction by exploiting coding region sequences of mRNAs.

MicroRNA are small non-coding RNAs that post-transcriptionally regulate the expression levels of messenger RNAs. MicroRNA regulation activity depends on the recognition of binding sites located on mRNA molecules. ComiR is a web tool realized to predict the targets of a set of microRNAs, starting from their expression profile. ComiR was trained with the information regarding binding sites in the 3'utr region, by using a reliable dataset containing the targets of endogenously expressed microRNA in D. melanogaster S2 cells. This dataset was obtained by comparing the results from two different experimental approaches, i.e., inhibition, and immunoprecipitation of the AGO1 protein--a component of the microRNA induced silencing complex.In this work, we tested whether including coding region binding sites in ComiR algorithm improves the performance of the tool in predicting microRNA targets. We focused the analysis on the D. melanogaster species and updated the ComiR underlying database with the currently available releases of mRNA and microRNA sequences. As a result, we find that ComiR algorithm trained with the information related to the coding regions is more efficient in predicting the microRNA targets, with respect to the algorithm trained with 3'utr information. On the other hand, we show that 3'utr based predictions can be seen as complementary to the coding region based predictions, which suggests that both predictions, from 3'utr and coding regions, should be considered in comprehensive analysis.Furthermore, we observed that the lists of targets obtained by analyzing data from one experimental approach only, that is, inhibition or immunoprecipitation of AGO1, are not reliable enough to test the performance of our microRNA target prediction algorithm. Further analysis will be conducted to investigate the effectiveness of the tool with data from other species, provided that validated datasets, as obtained from the comparison of RISC proteins inhibition and immunoprecipitation experiments, will be available for the same samples. Finally, we propose to upgrade the existing ComiR web-tool by including the coding region based trained model, available together with the 3'utr based one.

RIP-Chip analysis supports different roles for AGO2 and GW182 proteins in recruiting and processing microRNA targets.

BACKGROUND: MicroRNAs (miRNAs) are small non-coding RNA molecules mediating the translational repression and degradation of target mRNAs in the cell. Mature miRNAs are used as a template by the RNA-induced silencing complex (RISC) to recognize the complementary mRNAs to be regulated. To discern further RISC functions, we analyzed the activities of two RISC proteins, AGO2 and GW182, in the MCF-7 human breast cancer cell line. METHODS: We performed three RIP-Chip experiments using either anti-AGO2 or anti-GW182 antibodies and compiled a data set made up of the miRNA and mRNA expression profiles of three samples for each experiment. Specifically, we analyzed the input sample, the immunoprecipitated fraction and the unbound sample resulting from the RIP experiment. We used the expression profile of the input sample to compute several variables, using formulae capable of integrating the information on miRNA binding sites, both in the 3'UTR and coding regions, with miRNA and mRNA expression level profiles. We compared immunoprecipitated vs unbound samples to determine the enriched or underrepresented genes in the immunoprecipitated fractions, independently for AGO2 and GW182 related samples. RESULTS: For each of the two proteins, we trained and tested several support vector machine algorithms capable of distinguishing the enriched from the underrepresented genes that were experimentally detected. The most efficient algorithm for distinguishing the enriched genes in AGO2 immunoprecipitated samples was trained by using variables involving the number of binding sites in both the 3'UTR and coding region, integrated with the miRNA expression profile, as expected for miRNA targets. On the other hand, we found that the best variable for distinguishing the enriched genes in the GW182 immunoprecipitated samples was the length of the coding region. CONCLUSIONS: Due to the major role of GW182 in GW/P-bodies, our data suggests that the AGO2-GW182 RISC recruits genes based on miRNA binding sites in the 3'UTR and coding region, but only the longer mRNAs probably remain sequestered in GW/P-bodies, functioning as a repository for translationally silenced RNAs.